Spin precession in disordered systems: 
Anomalous relaxation due to heavy tailed field distributions 
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We investigate spin precession in the presence of randomly distributed field sources. Their fields 
hi reorient by thermally activated transitions and decrease as r _M with the distance from the spin 
probe. Based on analytical calculations and scaling arguments we find that the polarization decay 
of a spin ensemble exhibits a rich behavior characterized by stretched exponenentials and power 
laws with exponents depending on /i and the dimension d. The anomalous relaxation laws result 
from heavy tailed local field distributions and are verified by computer simulations. Implications 
for experiments are pointed out. 
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Many experimental probes rely on a precession of a 
spin S in an external field H, 



dS 
~dt 



S x H. 



(1) 



Examples are nuclear and electron magnetic resonance 
(NMR, ESR), muon spin relaxation (fiSR), /3NMR, and 
quantum optical measurements, where transitions in 
two-level systems can effectively be described by an equa- 
tion of type (|j) . In disordered systems the field H gener- 
ally exhibits both spatial and temporal fluctuations and 
the relaxation of an initially polarized spin ensemble is of 
interest. While traditionally this relaxation dynamics is 
studied for Gaussian stochastic processes H(t) [0J2), more 
complex stochastic processes became of interest recently 
(see e.g. ||). Here we will focus on systems, where the 
second moment (H 2 ) of the field distribution diverges. 
These situations occur, when the field H = ^ results 
from randomly distributed sources i in d dimensions with 
a spatial field dependence hi ~ r~ M , /i > d/2 (for dipolar 
fields, in particular, /i = 3). 

As an example of practical importance we focus on 
/iSR in disordered systems of single domain ferromag- 
netic particles In these systems the clusters per- 
form thermally activated transitions between certain easy 
magnetization directions with a rate v, which lead to 
fluctuations of the magnetic field at the muon site. We 
will show in this Letter that these fluctuations give rise 
to a rich anomalous relaxation behavior due to the fact 
that the random spatial distribution of the cluster mo- 
ments leads to Levy type local field distributions. De- 
pendent on how the reorientation rate v compares with 
the characteristic width W of the field distribution and 
dependent on the number of possible orientations of the 
cluster moments, we find very different relaxation sce- 
narios. The long time relaxation is given by either power 
laws or stretched exponentials, where the exponents de- 
pend on both fj, and d. The slow relaxations occur even 
in the absence of cluster interaction effects and in this 
respect should be contrasted to the relaxation found in 



spin glass systems || or related disordered systems |J 
with strongly interacting components. 

To be specific, we consider the following model. We 
place a spin S at the origin of a d-dimensional system 
that contains randomly oriented point-like clusters with 
number density n at random positions. A cluster with 
moment m and position r is assumed to induce a field 
contribution h = m/r^ at the probe site. Each moment 
m changes its orientation to a set of possible other ori- 
entations with the rate v. In particular we study two sit- 
uations: In the first case only the directions m and m 
are possible (uniaxial case), while in the second case there 
are four additional orientations perpendicular to m corre- 
sponding to a cubic symmetry (multiaxial case) . Initially 
the spin is polarized in the z-direction, S = (0, 0, 1). The 
task is to solve eq. (|l|) for a given cluster configuration 
and a certain realization of the cluster reorientation pro- 
cess and to average this solution over all possible realiza- 
tions. By finally averaging over all cluster configurations 
we obtain the spin polarization (S z (t)) at time t as mea- 
sured in experiment. In the following we will discuss the 
relaxation behavior for the generic situation [i>dj1 

We start out by focusing on the time regime t <C v~ , 
where the field H can be viewed to be static, and 
the solution of eq. @ reads S z {t) = {Hi/ H 2 ) + [1 - 
(H 2 /H 2 )]cos(Ht). By an exact calculation we obtain 
for the probability density ip(H) of the local field H 



tf(H) = 



(2) 



where ReL' a {u) denotes the real part of the derivative 
of the Levy stable law L a .o(u) = (2tt)~ 1 J dkexp(—iku — 
\u\ a ) to the index (a, 0) (see e.g. ||); the characteristic 
width W — Cw"tnn il l d is given by the field associated with 
the mean distance n~ 1//d of the clusters times a constant 
g. For large H, 4ttF 2 V'(H) - C^W^^H /W)- 1 - 1 '^ , 
implying that (H 2 ) does not exist. Averaging S z (t) over 
-0(H) eventually yields 
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FIG. 1. Spin polarization (S z (t)} as a function of ut in the slowly fluctuating case (v/W = 10 -3 ) for (a) multiaxial and 
(b) uniaxial cluster moments, and several {i and d. The symbols refer to the simulations and their assignment is the same in 
both figures. The dashed lines refer to the exact result (pf), while the solid lines are fits according to the long-time behaviors 
(^,^). The inset in (a) shows, on a semi-logarithmic scale, the exponential long-time relaxation of (S z (t)) vs. vt that is almost 
independent of d and [i (the solid line is drawn as a guide for the eye). The inset in (b) demonstrates the scaling (^|) for 4 
different radii n < n~ 1/d , n = 1.0 (+), 1.5 (x), 2.0 (*), and 2.5 (o) in the case fx = d = 3, n = 0.01. 



(S z (t)) 



■-(Wt) 



d/fi 



exp 



-(Wt) 



d/ii 



(3) 



For d — /i, i.e. in particular for dipolar fields in d = 3, 
one recovers the Lorentzian Kubo-Toyabe function |Q. 
As shown in Fig. [l] for different fi an °- d, the results from 
our simulations agree with eq. (||) for vt <C 1. Laws 
of type (Q) have been used in the literature to describe 
anomalous /iSR line-shapes with d/fj,^= 1,2 that neither 
follow a Lorentzian (d/fi=l) nor Gaussian (d//x = 2) be- 
havior (see e.g. [0). We note, however, that (||) is an 
exact result and should not be confused with an effec- 
tive "power Kubo-Toyabe function" [jll| that serves as a 
fitting function. 

In the dynamic regime t 3> v~ x we distinguish between 
the two cases of slowly or rapidly fluctuating cluster mo- 
ments, where v <C W or v 3> W, respectively. In both 
cases we employ scaling arguments to derive the typi- 
cal decay rates T of the spin polarisation. To tackle the 
problem of averaging over spatial cluster configurations, 
we consider subensembles of configurations that are spec- 
ified by fixing the distances of the clusters closest to the 
spin probe. This concept is motivated by the hierarchy 
implied by the Levy statistics, which for the field dis- 
tribution (||) means that the nth nearest cluster gives a 
contribution of order ?W d times smaller than the closest 
cluster (see e.g. fl2|). 

Let us begin with the case v -C W of slowly fluctuat- 
ing cluster moments, where for the relevant cluster con- 
figurations the field H has a magnitude H ^> v (other 
configurations have an exponentially small weight). In a 
time interval of order v~ x then, the spin precesses many 
periods around the local field, whereby S z (t) oscillates 
around a mean value S z (t). The changes of S z (t) av- 
eraged over many realizations of the cluster dynamics 
determine the decay of spin polarization. 



In the multiaxial case, significant changes of H, which 
occur in a time of order alter the axis of precession 
and S z (t) relaxes with a rate proportional to v. Hence 
we expect a simple exponential decay 



(S z (t)) ~ exp(-cst.vt) . 



(4) 



which is confirmed by our simulations shown in Fig. |l]a. 

The uniaxial case is more subtle. To see this, we de- 
compose the field H into the contribution hi = m/r^ 
from the nearest cluster at distance r\ and the contribu- 
tion Hi from the other clusters, H = hi 4- Hi. In the 
subensemble of all cluster configurations with given n, 
the variance of Hi is 



(Hl\r x ) = C H hl 



-d/fi 



(•5) 



For n ^>n _1 / d , hi/W <C 1, and Hi dominates over hi. 
Hence one encounters the same physical situation as in 
the multiaxial case. For small r\ <C however, hi 

is dominant, so that changes hi — > hi essentially revert 
the direction of precession and leave S z (t) unchanged. 

In this situation of small r x <C nT x l d the presence of 
the contribution Hi causes the axis of the field H (irre- 
spective of its direction) to wobble around the ±hi-axis 
with the rate v and an angular amplitude of order Hi /h±. 
The wobbling motion together with the much faster pre- 
cession leads to a diffusive type of motion of S z (t) with 
a diffusion rate T ~ [H\jh\) 2 v . 

To extract the asymptotic relaxation of the spin po- 
larization we consider the subensemble of all cluster con- 
figurations with fixed distances n and r-i of the near- 
est and second nearest cluster to the spin probe. In 
the configurations of this subensemble we can decom- 
pose Hi into h2 and H2, where /i2 = m/r2 and (-fffks) 
satisfies (H) with h\ replaced by h 2 - Accordingly, for 
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ri < r 2 < U\ ~ m 2 /r^ and T = T(r 1 ,r 2 ) oc 

(r 1 /r 2 ) 2f *is, while for r^^- 1 ^, £T?~ W d ^(m/r^) 2 ~ d l ^ 
and r(n,r 2 ) oc r^*(m/rJ) 2 - d /V. Writing 
(S z (t)\ri, r 2 ) ~ exp[— r(ri, r 2 )t] in the subensemble with 
given r\ and r 2 , we can average over the probability 
density 2 (r 2 |ri) = Sdnrj' 1 exp[— Vdn{r d — rf)} of r 2 



(ri < r 2 < co) to obtain |13| 



(S.WIri) ~ exp J V d nri-cst. Un^r^ut 



d/2ju 



(6) 



for i4 3> 1 (and r\ <Cn _1 / d ). We have verified this predic- 
tion for various /i and d by our simulations. One example 
(for fi = d — 3) is shown in the inset of Fig. [l]b. 

Final averaging over the probability density 0i(ri) = 
Scin-rf' 1 exp[— Vdnrf] of r\ yields 



(7) 



This slow power law decay is in marked contrast to the 
exponential decay in the multiaxial case and it is verified 
in Fig. 03 by our simulations. 

Next we discuss the case v 3> W of rapidly fluctuat- 
ing cluster moments. The field H in the relevant cluster 
configurations now has a magnitude H <C v and the spin 
rotates only by a small angle in a time interval of order 
v~ x . This means that the concept of a mean value S z (t) 
is not useful any longer, since the phase of the precession 
matters. Reorientations of hi are effective for the spin 
relaxation both in the presence of uniaxial and multiaxial 
cluster moments. 

The small angular changes of the spin lead again to a 
diffusive type of motion of S z (t). In time v the angular 
change is of order H/v and the corresponding diffusion 
rateT ~ (H/v) 2 v. Decomposing the field H = hi + Hi as 
before, and taking into account the dominant contribu- 
tions we thus find T = T(ri) (xv~ 1 m 2 /r 2 ^ for ri<n~ 1 / d 
and r(n) oc v- 1 W d ^{m/r>f) 2 - d ^ for n > n -^ d [cf. 
eq. (g)]. We then write {S z (t)\n) ~ exp[-r(n)t] for 
24 » 1 and r x > (m/u) 1 ^ (for n < (m/v) 1 ^, hi > v, 
i.e. one encounters a situation corresponding to the case 
of slowly fluctuating cluster moments). This exponential 
decay of (S z (t)\ri) is demonstrated in the inset of Fig. ^ 
for d=fi= 3 in the regime n >n~ 1 / d . By averaging over 
ri we finally obtain 



(S z (t)) ~exp 



—est. 



(v- x W 2 t 



d/2fi 



(8) 



To perform the average we have used a saddle point ap- 
proximation, where analogous comments apply as given 
in ]13] ], Figure |2| confirms both the scaling with {W 2 t/v) 
and the stretched exponential decay for the same \x and 
d values as in Fig. nl In the uniaxial case the stretched 
exponential decay (0) will, at long times, be masked by 
the much slower power law decay (p) that stems from the 
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FIG. 2. Spin polarization (S z {t)) as a function of W 2 t/v 
in the case of rapidly fluctuating cluster moments [u/W = 10 
(□), 50 (+), and 100 (x) for d = fx = 3, and u/W = 10 for 
the three other combinations of d and fj\. Data points refer 
to the simulations and the solid lines are fits according to 
eq. (^). The inset shows the exponential decay of (S z (t)\ri) 
and the scaling as discussed in the text for 4 different radii 
ri > n- 1/d , n = 6.5 (+), 7.0 (x), 7.5 (*), and 8.0 (o) in the 
case j« = d = 3, and n — 0.01 (the solid line is drawn as a 
guide for the eye). 

In summary we have shown that spin precession in 
the presence of randomly distributed and fluctuating field 
sources leads to an anomalous relaxation of an initially 
polarized spin probe, which is characterized by stretched 
exponentials [eqs. (||,|])] or power laws [eq. (0)]. The devi- 
ation from a simple exponential decay are caused by Levy 
type local field distributions [eq. (g)]. These render a 
treatment in terms of Gaussian processes impossible but 
allowed us to perform an analysis based on subensemblcs 
of cluster configurations that are defined with respect to 
the most dominant contributions to the local field, i.e. 
the field sources closest to the spin probe. 

It is important to stress that a simple mean field type 
description of the relaxation process would fail, as it was 
already pointed out by Uemura et al. p| in the case 
fj, = d — 3. In such a mean field description one might 
employ a "strong collision approximation" E3 , where 
the field H at the probe site is drawn anew from (||) 
with the rate v (thereby neglecting the fluctuations in 
the spatial cluster configurations) . By scaling arguments 
similar to those outlined above one can show that this ap- 
proach leads, for t 3> v , to an exponential relaxation 
(S z (t)) ~ exp(— r m ft) both in the cases of slowly and 
rapidly fluctuating cluster moments and irrespective of 
whether the clusters posses only one easy axis or more. 
For v <C W, one obtains r m f oc v, while for v 3> W, 
T nS Kv{Wlv)*ti> @. 

We restricted our treatment here to point clusters with 
unique moment m and neglected interactions between 
the moments. As long as the cluster sizes are much 
smaller than the mean distance nT 1 / d , crossover effects 
to a Debye like relaxation behavior typical for Gaussian 
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processes should be of minor importance. A broad dis- 
tribution of cluster sizes, however, may require a refined 
analysis in the dynamic regime (in the static regime the 
results remain unchanged except that m in the width W 
has to be replaced by its average value). To capture the 
dominant contributions to the local field and to take into 
account the variation in the jump frequencies (associated 
with changes in the anisotropy energy), it can be neces- 
sary to define the subensembles with respect to both the 
distance of the clusters nearest to the spin probe and the 
size of the clusters. Effects due to dispersion in the jump 
frequencies have been observed, for example, by /xSR in 
colossal magnetoresistive manganites [|l6j. Nevertheless, 
the basic scaling arguments presented in this work would 
still be applicable and an extension to systems of clusters 
with differing moments should be straightforward. 

Interactions between the cluster moments at high tem- 
peratures T can be accounted for by a temperature de- 
pendent width W — W(T) in (||) (for an approximate 
calculation in fi = d = 3, see |L7j). At low temperatures 
T by contrast, the cluster dynamics cannot be described 
any longer by a Poisson process with rate v (for dipolar 
systems in d — 2,3 this occurs for T < 0.5 m 2 n d ^ 3 , see 
jl8],[l9|] ) . In this low-temperature regime the problem be- 
comes more difficult and the relaxation laws (0,§) may 
no longer hold true. A non-Poissonian cluster dynamics 
has recently been encountered in a spin glass also ||. 

Having mentioned these limits of our findings, we hope 
that our work will stimulate further research on the chal- 
lenging problem of spin precession in disordered systems. 
Our scaling methods should give deeper insight into the 
spin relaxation in disordered systems and may be ex- 
tended to describe /iSR (or /3NMR) in other complex sys- 
tems, as e.g. spin glasses, structural glasses, amorphous 
magnets or disordered superconductors. 
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